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The behaviour of the master field in "induced QCD" near the edge of its support is studied. 
An extended scaling domain, where the shape of the master field is a universal function, is 
found. This function is determined explicitly for the case of dimensions, close to one, and 
the D — 1-expansion is constructed. The problem of the meson spectrum corresponding 
to this solution is analyzed. As a byproduct of these calculations, a new, explicit equation 
for the meson spectrum in "induced QCD" with a general potential is derived. 
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1. Introduction. 

The "induced QCD" is a lattice gauge theory of a scalar matrix field defined by 
the functional integral [[TJ 

Z= yp$(x)Pt/ M (x)exp|-iV^[trt/($(x)) 

-tr + fia)U^(x) 

In this formula a: marks the sites of the D-dimensional cubic lattice with lattice spacing 
a, are N x N Hermitian matrices and U^(x) is the gauge field (/i referring to the 

directions of lattice links, fx = 1, . . . , D). 

The field interacts with itself through the potential U(<&) and has the obvious 

kinetic term tr [$(x)L r M (x)$(x + p^U^x)]. This term is the only place where the gauge 
field enters this theory. Indeed, no special action for the gauge field is included in the 
model. 

Since U^(x) is incorporated in the action in such a simple way, the "induced QCD" 
is, in principle, exactly solvable in the large iV limit. In this respect ( |1 . 1| ) is one of the 
very few matrix theories which admit an exact solution not only in D < 1, but also in 
D > 1 dimensions. Needless to say, it would be extremely interesting to investigate such 
a solution, especially in the D > 1 case. In particular, the scaling properties of this model 
near its phase transition point must exibit universal features, characteristic of other D > 1 
models. 

In the large N limit the eigenvalue density p((p) of the field obeys classical 

equations of motion. That is to say, p(<p) plays the role of the master field in "induced 
QCD". Its classical dynamics can be formulated in two ways. One way is to use the 
Schwinger-Dyson equations for the functional integral ( |1 . 1] ) , obtaining a nonlinear singular 
integral equation for p((f)) 0. Quite remarkably, for the case of quadratic U(<&) this 
equation can be solved exactly ||. On the other hand, the same density p(<p) is related to 
solutions of a certain quasilinear partial differential equation (the Hopf equation) j| . These 
two formulations may appear very different, but in fact they are equivalent. Under some 
general assumptions about the potential U(<&) both of them reduce to the same functional 
equation !,§,§. 

In the vicinity of critical points, where p((p) vanishes, the small fluctuations of the 
eigenvalue density become significant, leading to critical behaviour. The location of the 
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critical point depends on parameters of the potential, but the critical indices are, presum- 
ably, universal. Remarkably, in the critical domain the master field equation simplifies, 
which allows one to determine the critical indices without constructing the complete solu- 
tion of the whole master field equation. 

The critical point can be positioned either inside of the support of p(<p) or at the edge 
of this support. The critical indices, corresponding to these two cases, are different. If 
p(<p) vanishes inside of its support (for an even potential this would occur at = 0) the 
eigenvalue density behaves at the critical point as 

p((p) OC |</>| 1+7 , COS7T7 



3D - 2 

This case is rather sophisticated. In particular, an important problem of the mass spectrum 
corresponding to this critical point is still unsolved. 

The other, simpler situation was considered recently by Boulatov 0. He assumed 
that the master field vanishes at the endpoint of its support according to 

p(4>) oc (0-a) 1+ ^(0-a). 

He found two possible solutions for 7: 

D 

cos 717 = D or cos 7r7 
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However, it is still unclear how the Boulatov's master field behaves away from the 
endpoint. The scalar field potential giving rise to this solution is not known either. Finally, 
although, by construction, this master field corresponds to an extremum of the action of 
the "induced QCD" , it may maximize, rather than minimize, the action. 

In this paper we will take a closer look at Boulatov's solution. In section 2 we de- 
termine the shape of the master field away from the endpoint. We find that there is a 
whole domain where this field is described by certain real universal functions, which can 
be represented in terms of a power series. In section 3 we consider the case of dimensions, 
close to 1, and construct the D — 1 -expansion of the master field. This allows us to sum 
the series and evaluate these universal functions explicitly. 

Finally, we address the problem of spectrum in "induced QCD". In section 4 we use 
the connection between the "induced QCD" and the Hopf boundary problem to obtain a 
new, simpler form of the wave equation, describing the meson spectrum of the theory. This 
derivation applies to "induced QCD" with any potential and is not in any way restricted to 
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Boulatov's solution, the spectrum of which we shall determine in section 5. Our calculations 
show that in this case the meson spectrum contains tachyons. This implies that such 
solution does not correspond to a local minimum of the action. Nevertheless, the master 
field we are investigating exibits some general features which are likely to be present in 
other, healthier solutions. 



2. The Series Expansion of the Master Field. 

The exact solution of the large- N "induced QCD" is based on the saddle point method 
0, 0. First, one integrates out U^x) in using the Itzykson-Zuber integral || 

/(<&, *) = J VUexp[Nti$U*rf) = det ^^ j) ^ (2- 1 ) 

where 

4>i and ijjj being the eigenvalues of the N x N matrices $ and One is left then with an 
effective theory of the "eigenvalue fields" <j>i(x) which has the actionB 



Seff [<£(*)] =Nj2 U Mx)) ~ 

x,i 

£ln[l(0(aO,#c + H) 



2.2 



+ 



In the large iV limit the remaining functional integral over (pi(x) is dominated by the saddle 
point of S e s[(f)(x)], that is, by x-independent <fii(x) = <pi satisfying 

d -S eS [cf ) ] =0. (2.3) 



Since the effective action S e R is invariant with respect to any permutation of <f>^i = 
1, . . . , N, the saddle point can be fully described by specifying the density of eigenvalues 
(pi on the real line, p(<f>). This density plays the role of the "master field" in this theory. 
To solve the "induced QCD" means to find p(<fi) for the given interaction potential U(<&). 
This, in turn, allows one to evaluate other physical observables of the theory. 



1 The second term in S e s is due to the fact that the integration measure T>&(x) after the 
change of variables to 4>i(x) reduces to A 2 ((f>(x)) FJ^ T>(j>i(x). 
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As a consequence of the saddle point condition (pT3|), p(<f>) obeys a certain nonlinear 
singular integral equation 0. However, for a generic potential U ($) it is possible to derive 
a Sunctional constraint fixing „ W §■ One introduces the two functions^ 

G ± ( X ) = ± ) U'( X ) + ^pJ^±iM X )- (2.4) 

Then the saddle point of the "induced QCD" is determined by the equationsi 

G + (G-(x)) = x, 
G-(G+(x)) = x. 



(2.5) 



That is to say, G+(x), as a function of x, is inverse with respect to G-(x). 

Generally speaking, the master field p(x) is not equal to zero only inside a finite 
interval (referred to as the support of p(x)), and vanishes at its endpoints. However, it is 
the behaviour of p(x) in the vicinity of endpoints that determines the universal properties 
of the theory. For example, if p{x) ~ (a — x) 7+1 near the endpoint a, then the exponent 
7 is universal. 

The general solution of equations ( |2.5| ) is unknown. Therefore, there is no direct way 
to classify all possible values of 7 that can be achieved in "induced QCD" . Nevertheless, 
in the vicinity of a point where p(x) vanishes, the equations should simplify. Boulatov || 
has observed that (|2.5| ) admits a very simple solution, G+(x) = G-(x) = — x. Clearly, 
such solution by itself does not describe the "induced QCD" with any potential, since it 
contains no imaginary part at all. To correct this drawback, Boulatov perturbed G± with 
a powerlike term, 

G+(x) = -x + a(-x) 1+1 + . . . . (2.6) 
If 7 is not an integer, then for x > 

itp(x) — ImG+(x) = ax 1+7 sin(7T7) + . . . ^ 0. 

Boulatov found that this ansatz is consistent with ( |2.4[ ) only if 7 assumes a certain value, 

cos7T7 = D or cos7T7=— . (2.7) 

' 1 2D-1 K ' 



2 This formula defines G±(x) for real x. For complex x, G±(x) are defined by means of analytic 
continuation from the real axis. 

3 It is not known, however, whether or not some nonanalytic solutions of induced QCD can be 



described by this constraint. 
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In other words, the behaviour of the master field near the endpoint of its support is 
completely determined by the dimension D. Notice that, except for some special cases, 7 
is not a rational number. 

Since the constraint (|2.5|) is nonlinear, any perturbation like (— x) 1+7 will inevitably 
generate perturbations of higher orders. One ends up with a whole series for G+(x), 

00 

G+(x) = -x + Y,9k(~xy k+1 . (2.8) 
fc=i 

If 7 were rational, 7/c + 1 could become an integer for a sufficiently large k. In this case 
the k-th order term should be interpreted as induced by the (polynomial) potential U (x) 
in (|2.4| ). However, if 7 is not rational, such phenomenon can never occur, and all of the 
coefficients can be determined dynamically, from the master field equations ( |2.5| ) alone. 
This suggests that the whole infinite series ( |2.8| ), not only its first term (— x) 1+7 , has some 
universal properties. Let us show that this is indeed the case. 
To begin with, let us note that the function 

«*> = wn-nVM + / — (2 ' 9) 

2(D-1) J x-y 

is an analytic function of x, which has a cut along the support of p(x). Since on the cut 

f(x ± iO) = 1 U'(x) +vf P -^- T ^7rp(x), 
2(D-1) J x-y 

it follows from (|2]|) that 

G+(x) = -x+ 2D ~ 1 J(x) - —f(x), 
+ K ' 2D Jy ' 2D J K h , 10 x 

2D - 1 1 - 1 ' 

G- (x) = —x H — — fix) — fix). 

\ > 2D y ' 2D K ' 

Boulatov's ansatz ( |2.6| ) essentially means that f(x) has a branch cut along the real 
axis from to +00. Consequently, it can be expanded as 



00 



f( x ) = x J2f k (-xy k . (2.11) 



k=i 



This expansion makes sence by itself for real negative x. To define f{x) for real positive 
x, we need to perform the analytic continuation, with the result 



00 

tkl"x 

k=l 



f(x)=xJ2f^ k - lk 
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where 



Then 



z = e l7r7 . 



■ Z ~ k x lk . 



f(x) = x ^2f k i 
k=l 

Substituting this into ( |2.10|) and introducing the notation 

_ _ 2D -I _ k 

k ~ 2D 2D Z ' 

B 1 2D -K jk 
k ~ 2D 2D 1 

we easily derive the following expansions: 

—for real positive x on the upper edge of the cut, Imi = -HO: 



(2.12) 



G+(x) = -xfl + f^fkAkX^ (2.13a) 
^ fe=i ' 

G- [x] = -x(l + jr f k B k z- k x^\ (2.136) 
^ fe=i ' 

— for real negative x (there is no cut at x < 0): 

G+(x) = -x(l + J2fkAkZ k (-xy k \ (2.14a) 
^ fc=l ' 

G-{x) = -x(l + jTf k B k (- X y k \ (2.146) 
^ k=i ' 

Our goal is to determine f k . To this end, we must impose the constraint ( |2.5| ). This 
can be done in two different inequivalent ways. The first way is to impose G+ (G-(x)) = x 
on the upper edge of the cut, where Imx = -HO. 

In view of fl2.13a|) and ( 2.14£| ) this requirement can be written down as 



OO y OO OO \ 

(l + ^f k A k x^U + ^f p B p x^{l + ^f q A q x^y P ) =1. (2.15) 

fc=l ^ p=l 9=1 ' 

Then one can check that the other constraint, G-(G+(x)) = x, is fulfilled on the lower 
edge of the cut, where Imx = — iO. In addition, both of the constraints ( |2.5| ) are satisfied 
on the real axis outside of the cut, at x < (see fig. 1). 
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Fig. 1: (left) The functional composition G+ \G-(x)) takes a real point x > to 
itself. More precisely, G+ (G-(x)j maps the upper edge of the branch cut onto the 
lower edge, (right) Both G+{G-(x)j = x and G- (G+ (x) J = x are satisfied on 
the negative part of the real axis. 



The second way to impose ( p75| ) is to demand G + [G-(x)) = x on the lower edge of 
the cut. This will lead to an equation, different from ( |2.15| ). As the treatment of these 
two cases is completely analogous, we will now concentrate on the first one. After that, 
we will briefly describe the second case. 

Equation ( |2.15| ) allows us to find the coefficients fk recursively, one after another. To 
do this, we expand the right hand side of ( |2.15|) in a power series of a formal parameter 
t = x 1 : 

1 + (A x + B 1 )f 1 t + ((A 2 + B 2 )f 2 + (1 + 7 )Aj5 1 / 1 2 )t 2 
+ [(A 3 + B 3 )f 3 + ((1 + 1 )A 2 B 1 + (1 + 2 1 )A 1 B 2 )f 1 f 2 

Since fi ^ 0, we obtain the consistency condition 

A 1 +B 1 = 

as well as 

h 



t 6 + ... = i. 



"(1 + 7)-^/?, 



A 2 + B 2 



(2.16) 



[2.17a) 
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A = -( 1 + ^PTTTOT^( (1 + 7)( ^ B2 + B ^ ) 

^ A t n A \ i ^3 



+ ^(^2-^))^, (2.176) 



Using (|2.12|) and (|2.16| ), we see that the consistency condition translates into 



1 2D , 



which, in view of z = e 171 " 7 , implies the Boulatov's result (|2.7| ). Since 7 has to be real, this 
solution makes sense only for D > 1 or D < 1/3. 

We can simplify ( |2.17| ), using ( |2.12| ) to express Ak and Bk in terms of z, and eliminating 
D by virtue of ( gT§ ): 

f 1 + 7 f2 n , 2^ -1 2 

f (1 + 7K2 + 37) 3 (2D-1) 2 3 ^ 19 ^ 

/3 = 2(*-s + + 2 + z + Z 2) h - (1 + 7)(2 + 37) 4D(4D - 1) A ' 

and so on. Note the remarkable feature of these expressions: they are all real. Indeed, 
since z is a complex number of absolute value 1, the denominators of fk are invariant 
under complex conjugation. It is not at all obvious from ( |2.15| ) that this should be the 
case. Indeed, the reality of fk is closely connected to the fact that outside of the cut the 
first of the two equations (|2.5|) , G+ [G- (x)) = x, implies the second, G- (G+(x)) = x, and 
vice versa. Let us now prove that all of the fk are given by real numbers. 

We will do this in two steps. The statements below are easy to infer by inspection of 

When viewed as functions of {A±, A 2 , . . .}, {Bi, £? 2 , . . .}, the functions fk ({Ai}, {Bj}) 
are symmetric with respect to the interchange of A and B: 

fk{{A l },{B J })=fk{{B l },{A J }), 



provided that the consistency condition ( ]2.1b]) is satisfied. 



At first sight, this might appear false. Indeed, the expression for ^3 in (|2.176|) is not 



at all symmetric under the interchange of A and B. However, its asymmetric part equals 



and vanishes due to the consistency condition (|2.16| ). 

To prove the symmetry of fk, we introduce an auxiliary variable 



u = x(l + J2fkA k xt k ). (2.20) 

k=i 

Then it is a consequence of ( |2.15| ) that 

oo 

x = u (l + J2fp B P ulP )- ( 2 - 21 ) 

p=l 

Substituting ( |2.21| ) back into ( |2.2C1| ), we get 

OO y OO OO \ 

{1 + E f P B P ulP ) ( 1 + E fi A « uiq i 1 + E fiBi^y ) ( 2 - 22 ) 

p=l ^ q = l 1 = 1 ' 



u = u\ 



which is the same equation as ( [2.15| ), but with all replaced by and vice versa. Since 



the fk can be recursively determined from either of the equations (|2.15[ ) or ( |2.22| ) and since 
fk are unique, they have to be symmetric in A and B. 

If we assign to Aj and Bj a formal degree of j, then fk has degree zero for any k. 

That is to say, if we scale 

Af ] = X l A l: 

A) ( 2 - 23 ) 
B\ X) = A l S l5 

then 

fk{{A^ ) }ABf ) })=fk{{A l }AB J }). 

This statement follows immediately from the observation that such rescaling amounts 
to a redefinition of expansion parameter x in ( |2.15| ), x 7 — > Ax 7 , which does not change the 
equations for /. 

Now we are able to prove that fk is real for any k. Indeed, we can eliminate D in 
( |2.12|) in terms of z, using ( |2.18|) . This gives 



B*W = 1-^(1 + ^'). 

Under complex conjugation, 

Al( Z ) = Z - k Bk(z), 
Bl{ Z )= Z - k A k {z). 



(2.24) 



Therefore, 

r k {{Ai}AB j })=fk{{A*},{B*})=fk{{z- i B i },{z-iA j }). 
Setting A = z~ x in ( ggg) , 

/*({*-'£,}, {ar'^}) = /*({Si},{i4i}) 
and, by the symmetry proprety, 

/ fc ({5 i },{A i }) = / fc ({A i },{5 J }). 

Hence 

/fe({^} ! {^})=/ fc ({A i },{B j }), 

so that all of the are indeed real. 

One can see from ( |2.19| ) that fk/fi are universal numbers, which depend only on 



the dimension D. The coefficient fx (which can be arbitrary) enters all formulas in the 
combination fix 1 , thus setting the scale for x. We conclude that there is a whole scaling 
domain, where fix 1 ~ 1. The (universal) behaviour of the master field in this domain is 
more complicated than the powerlike shape of the singularity at x — > 0. In the next section 
we will investigate this behaviour in more detail. 

Finally, let us note that the condition on 7, fl2.18[ ), admits not only a positive, but 
also a negative solution, 7 = —(1/tt) arccos _D/(2.D — 1). In this case ( |2.6| ) is not a good 
approximation at all. However, if one takes into account the whole series for f{x 1 ), it is 
possible to derive the behaviour of p(x) as x — > 0. In the next section we will see that in 
this case the master field develops a logarithmic singularity. 



3. The D — 1-expansion. 

One might wonder how the universal function / behaves at finite values of its argu- 
ment. Obviously, the first few terms of a Taylor series do not answer this question. One 
needs to evaluate all coefficients of the series and perform the summation. 

This can be done explicitly when 7^0. Indeed, it is easy to see from (|2.19|) that the 
fk(l = 0) are not singular. We will see below that in fact they describe the first order in 
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the D — 1-expansion around D = 1. To calculate the numbers f k = fk{l = 0), we return 
to the general formula ( |2.15| ), with t 



x 1 : 



OO s OO OO \ 

fc=l ^ p=l (J=l ' 

Let us take the limit of this equation as 7 — > 0, while keeping t fixed. Using ( |2.24|) we 
derive 

^-2 2 

A k (z) =iTT 7 k~ — -- + 0( 7 3 ), 

2 (3-2) 
S fe (z) = -m 7 /c + ttV (fc 2 --)+ 0( 7 3 )- 



Expanding flOl ) up to C(7 3 ) and keeping in mind that Ak(z) ~ B k (z) ~ 0(7), we get 

00 00 

+ Yl fkf P A k B p t k +P + O('y 3 ) = 0. 

k=l k,p=l 

By virtue of ( |3.2j ) this means that the /fe(7 = 0) satisfy the identity 

00 00 
Y,(k 2 - l)f k t k + kpf k f p t k+ ^ = 

k=l fc,p=l 

which translates into an ordinary differential equation for 

OO 

f(t) = J2fkt k , 

k=l 

Introducing a new variable £ = log£ and denoting / = df /d£, we obtain 

f + f 2 ~f = 0- (3-3) 



The solution of this equation is given by 

/(*) 



logt = / \ ^ du. (3.4) 

' Xl 2u + e~ 2u - 1 K J 



Taking into account that f(t) ~ fit as t — > 0, we finally obtain 

/(*) 



l ° g W) = / { V 2, + e-^-l " ^ V- (3 - 5) 
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oo 

E 

k=l 

^ (3.6) 



In the same approximation the master field equals 

Trp(x) =ImG + (x + i0) = -x^2f k t k lmA h 

k=i 

oo 

= - 7T7X kf k t k = -Tvjxf(t 
k=l 

= - Ti^x^fix-i) - i(l-e- 2 /(- 7 )). 

Note that we have fixed x 1 to be finite, although 7 — > 0. This means that in the scaling 
domain, where our expansion is valid, |logx| ~ 0(1/7). 



Since p(x) has to be positive, ( |3.6|) implies that f"f < 0. Therefore, there are two 
possibilities: either 7 < and / > 0, or 7 > and / < 0. In the second case the terms in 
the series expansion of /(£) get smaller and smaller, as x —>■ 0, and Boulatov's result gives 
a good approximation. This picture fails, however, in the first case, when 7 < 0. Then, 



as x — ► 0, we have to use the t = x 7 — ► 00 asymptotics for fit). From (|3.4| ) we derive, 
remembering that / > 0, 

/(*) 



so that 



and 



J V u 



f{t) ~ ^log 2 t, 



7 2 

pix) ~ — — xlogx, x — > 0. 



We see that in the scaling regime with 7 < the master field develops a logarithmic 
singularity at the endpoint of its support. 
Since as D — > 1, 

7 = ±- V / 2(D-l), 

TV 



( p.6|) represents the leading order of the yD — 1-expansion around D = 1. Note that the 
expansion parameter is \/D — 1, rather than D — l. This is related to the fact that at D — 1 
the master field vanishes identically. The perturbation, therefore, is singular in D — 1. 

Higher orders of expansion can be obtained in a similar way. To evaluate the first 
correction, we set f k = fkil = 0) + "if^ + • • • an d expand ( |2.15|) to C(7 4 ). After some 
algebra one derives that the function 



(i) + fe 

fc=i 
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satisfies an ordinary differential equation: 



/(i) + 2//<i> - /W 



where = df^/d^, £ = logt, and / is the leading order result given by (|3.4|). 



In fact, this D — 1-expansion gives a quite accurate approximation. We have calculated 
the function f(t) numerically for D = 2 by the following extrapolation technique. We used 
( |375| ) as a guide, and modified it as follows 



The function g{u) with (7(0) = 1 was computed as a Taylor series from exact equation at 
D^l. Then, the Pade approximation was used. We checked the precision by comparing 
diagonal with nondiagonal approximants. The 15/15 and 14/16 approximant practically 
coinsided for —2 < u < 2. This is where we numerically computed the corresponding 
integral, with 15/15 approximant for g(u). 

The result, along with the plot of the function / for 7 = is presented in the figure. 
We checked numerically that the difference remains small even at D = 3. 




(3.7) 







fit) 



2r 



1.5- 



-1 . 5 




-2 L 



1 . 5 



t 



Fig. 2: Plots of the function f(t) for D -> 1 (7 = 0) and for D = 2. 
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Analogous calculations can be carried out for the second solution, which is obtained 
by imposing the constraint G+(G-(x)) = x on the lower edge of the branch cut (see the 
discussion after fl2.15|) ). This leads to another solution, valid at \D\ < 1, 

z + - = 2D, (3.8) 

z 

It is possible to construct 1 — .D-expansion (in contrast with D — 1-expansion we discussed 



above), with 7 = ±a/2(1 — D)/n, np(x) = — 7T7/(x 7 ), and 

/(*) 

'°g'=/V ^-L-A ' 3 ' 9 » 

In this case, if 7 < 0, f(t) is defined only for fit < 1.27635, meaning that the ansatz ( |2.11j ) 
is internally inconsistent. The only self-consistent possibility in this case is 7 > 0, / < 0. 
In fact, the universal function defined by (|3T9|) is just a continuation of the function / 



in (|3.5| ) to t < 0. Indeed, using (|3.5|) and (|3.9|) , and keeping in mind that f(t) ~ f(t) ~ fit 
as t — > 0, it is easy to check that 

ht) = -n-t). 

We see that the case D > 1 is described by the function f(t) with t > 0, while its analytic 
continuation to t < will give the solution for D < 1. 

If any of these solutions is to describe the "induced QCD", it has to minimize the 



action in (|1.1|) . However, the master field equations only guarantee that the solutions 
above are the extremal points of the action. To find the type of extremum corresponding 
to these solutions one has to consider the quadratic variation of the action functional and 
determine the normal modes and eigenfrequencies of small oscillations around the master 
field. 



4. General Theory of the Spectrum in "Induced QCD". 

In this section we will obtain an equation describing the normal modes of small fluc- 
tuations around the master field in the "induced QCD". An integral equation answering 
this question has been derived before 0. However, we will derive a new, simpler, equation 
in a different way. The result will be applicable to the "induced QCD" with any potential 
U ((/)). We will apply it to the particular case of the Boulatov's solution in the next section. 
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As one can see from (|2.2|) , the effective action of "induced QCD" in the large iV limit 



is 



^eff=iV 2 ^| J p x {<f,)U{<l>)d<f,-^2F[p x {<f,),p x+lic 




rf0lrf02Pa;(0l)p a; (02)ln|01 - </>2 



(4.1) 



where p x {4>) ls the eigenvalue density at the lattice site x, and F[p((j>),cr(ip)] is the large 
N asymptotics of the Itzykson-Zuber integral: 

F[p(cf>),a^)] = lirn^ ^ln/[*,tf] 

p((f>) and cr(ifj) being the eigenvalue distributions of matrices $ and \&. 

We would like to see whether or not the small perturbations of the master field are 
unstable. To this end, let us derive the equation of motion for these small perturbations 
(we shall refer to it as the wave equation). 

Writing 

Px(<f>) = Po(<t>) + 5p x {4>) = Po{<f>) +^e im 8p$(4>) 



where po(<f>) is the translationally invariant master field, determined from (2^3), we find 
that the part of the action, linear in 5pp{4>), vanishes due to the master field equation 
(p73|), while the part quadratic in the fluctuations Sp^((f>) is 



S'jSeff = - ^ J #i# 2 <5Pp(0l)-^p(</>l,^2)<5P-p(</>2) 



where 



K 



pWl: 



b 2 ) = In |0i -0 2 | +D 



S 2 F[p,a] 



5p(4> 1 )5p((p2 



p(x)=cr(x)=p (x) 



2 _ 6 2 F[p,a] 



5cr(0i)5p(0 2 y 
D 



p(x)=a(x)=po (x) 



(4.2) 



(4.3) 



^ 2 (p) = cosa ?V 
p=i 

Note that since p(x) = cr(x), the kernel Kp is symmetric in <pi, <p2- 

Since J p((f))d<p always equals 1, the perturbations Spp(<p) must satisfy the normaliza- 
tion condition 

8pp((j>)d(f) = 0. 
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Therefore, we can introduce the functions Vyf^) such that 

Sppi<i>) = -^^M » ^(-oo) = Tpp(+°°) = °- 

The equations of motion for ipp(4>) can be obtained by varying with respect to Vy(<^) : 

= 2/^,^)^ = 0. (4.4) 

This equation has solutions only for special discrete values of O^, which form the spectrum 
of "induced QCD". 

To proceed, we will need to evaluate the second derivatives of the Itzykson-Zuber 
integral, which are present in the kernel ( |4.3|) . Although /($, is in general given by a 
simple explicit formula ( |2.1| ), the function F[p, a] has a very complicated structure. 

Let us first outline the calculation of the first functional derivatives of the Itzykson- 
Zuber integral. It is convenient to introduce the two functions v(<p) and u(<p), defined 
by 

d 8F[p{4>)M^)] =(j) + v ^ ) _ v MW 



d<t> 6p(</>) 



dtp Sa(tp) J ip — ip' 

Then v (<p) and u(tp) can be determined as follows J3J. One introduces an auxiliary complex 
function (p(x, t) and sets up the boundary problem for the so-called Hopf equation: 

d<p(x,t) d(p(x,t) 

lmip(x,t = 0) =7rp(x), ( 46 ) 
lm(p(x,t = 1) = na(x). 

Then it is possible to prove that 

v (x) = +Re (f(x, t — 0), 
u(x) = —Re(p(x,t = 1). 

Let us emphasize that that we are considering the case when the two eigenvalue distribu- 
tions entering the Itzykson-Zuber integral are not necessarily the same. 
In terms of v (x) and u(x) the kernel of the wave equation becomes 

= (i _ D) i + D pM + (4 . 7) 
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To find the functional derivatives 8v((f>i) / 'Sp((f)2) and 8u{<f>\) / '8 p((f>2) , which correspond to 
the second variations of the Itzykson-Zuber integral, we have to see how the solution of 
the Hopf equation (03) changes when we vary the boundary conditions. This will allow 
us to derive two integral equations, constraining 8v(<pi)/8p((f)2) and 8u((f>i) / 8 p(4>2) • The 
easiest way to do this is to use the fact that the Hopf equation has infinitely many integrals 
of motion. We can choose 

H &) = ~ 7rt fq + i) Im / dx(x-t<p(x,t)) q+ \ (4.8) 

It is a consequence of ( fOj) that dH q /dt = 0. Hence H q (t = 1) = H q (t = 0), which, in 
view of the boundary conditions for (p entails 



7T{q 



Considering a small variation of p(x) and setting Sa(x) = (that is, varying one of the 
distributions in the Itzykson-Zuber integral while keeping the other fixed), we get 

Im J dxSu(x) (x + u(x) — Z7ra(x)) 9 = J dxx q 8p(x). 

Summing these equations for q = 0, 1, . . . with weight z~ q ~ 1 , we obtain, provided that z 
is outside of the support of p(x): 

-hm I dx ; 5 M = / dx 5 ^. 



J z — (x + u(x)) + ina(x) J : - ■'' 



Now we can put a(x) = p(x) = po(x), u(x) = f(x), and, denoting R(x) = x + u(x) = 
x + v (x), we obtain the first integral equation, which determines 8u(x)/8p(y): 

3 fdx( 6u{x) ) np{x) - 3 - (49) 

dy J V Sp{y) ) (z- R(x)f + ix 2 p 2 {x) dy z - y 

To obtain an equation for 8v(x)/8p(y), we use another set of conservation laws 
Iq ® = ?r(l -t)(q + l) Im / dx ( x + ( 1 ~ t ^( x ^))' 



Since now 



^ = l)=/<«*) = 0, 
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the same procedure yields 

8I q (t = 0) = — Im J dx(5v(x) + in5p(x)) [x + v(x) + inp(x)) q = 0, 



so that 



hn / dx Sv(x) + inSp(x) = Q 
z — R(x) — iixp(x) 



and 







dx 



Sv(x) \ 



7rp(x) 







z-R(y) 



dy J " \Sp(y)J (z-R(xj) +tv 2 p 2 (x) dy (z- R(y)) +Tt 2 p 2 (y) 
Notice that the integration weight in ( |4.9| ) and ( |4.10| ) is the same, 

7Tp(x) 



(4.10) 



(z-R(x)) +n 2 p 2 ( 



X) 



So we can integrate the wave equation (|4.4|) over x with this weight to eliminate 
Sv(<pi)/8p((f)2) and 8u(<pi)/8p(<p2)- Keeping in mind (|4.7|), we get the final form of the 
wave equation^ 



z-y 



-D 



z - R(y) 



(z-R(y)Y + 7r 2 p 2 (y) 

dx ivp(x) 



+(D-1)V 



x — 



y (z-R(x)) +iv 2 p 2 (x) 



(4.11) 



8pp(y)dy = 0. 



The functions R(x) and p(x) are known from the solution of the master field equation. In 
fact, using ([4.5|), ( [Of ) and ( p.4|) , one can deduce that 

G±(x) = R(x) ± iirp(x). 

The equation we derived here is very explicit. It contains no auxiliary functions and 
can be written down at once for any given solution of the master field equation. For 
example, it can be used to derive the spectrum of the "induced QCD" with quadratic 
potential 0, recovering the result of Aoki and Gocksch 0. We will use this equation to 
analyze the Boulatov's spectrum. 



We remind that, by construction, z is positioned outside of the support of p(x). 
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5. The Spectrum of the Solution with the Endpoint Singularity. 

The functions R(x) and irp(x) in the Boulatov's case equal (see ( |2.6| )): 



R(x) = —x — a cos 7T7 x 1+1 + 



ixp(x) — asin7r7x 1+7 + 



As x — > 0, we can approximate the weight in (|4.11| ) 

np(x) 



R(x)) +7V 2 p 2 { 



5(z - R(x)) ~ 8{z + 



x) 



X 



In this approximation the principal value integral, present in the wave equation, can be 
estimated as follows 



V 



dx 



np(x) 



X ~V (z - R(x)) + n 2 p 2 (x) 
With the same accuracy, 



r j J^s(z + x) 
x-y 



1 



z + y 



z - R(y) 



V- 



V- 



(z-R(y)) 2 + n 2 p 2 (y) z ~ R (v) - + H 
and we see that Qlll|) in the leading order becomes 



> + <p. 



5pp(y)dy = 0. 



(5.1) 



z-y z + y_ 

This equation admits powerlike solutions, 5p^(y) = y a . To find Vt 2 {p), we use the formulas! 

y a dy , .„ 

= — 7rcot 7ra (— z) , 

(5.2) 



V 



z + y 

y a dy tv 



This implies 

At small lattice spacings 



z — y sin ira 



Q 2 (p) = —cos na. 



D 1 

il 2 {p) = cos ap^ = D- -a 2 P 2 . 



We are keeping in mind that z < 0. 



19 



We see that, if D > 1, whatever the index a is, the particles of this theory have P 2 = 
— m 2 > 0. Hence, m 2 < 0, that is, these particles are tachyons. 

To study the D < 1 case we have to determine the index a. We shall argue that 
a = 1 + 717, with n an integer. 

Indeed, the wave equation ([4.4j ) means that, as long as po(<fi) solves the master field 
equation, so does po(<fi) + Sp x (4>). Since po(4>) has an expansion in powers of </> 7fc , it is 
natural to look for small, coordinate-dependent, perturbations of the coefficients of this 
expansion. This suggests that the n-th normal mode of the of the wave equation arises if 
we perturb the n-th coefficient in, say, (|2.8|). Such mode will be given by 



5 {n) P P < 



x 



x 



l+7n 



\p) x 



k=0 



{x 1+1 x 1+2,1 



The functions 5^ Pp(x) can be obtained by linear combinations from the infinite set 
.}. A very important property of this set is that it remains invariant under 
the action of the wave operator dK/dx. 

Let us represent the action of the integral operator in ( |4 . 1 1| ) by an infinite-dimensional 
matrix M: 



7rp(x)dx 



dy 



dK(x,y) 
dx 



(Z-R( X )) +7T 2 p 2 (x) 

The matrix M has an upper-triangular form 



/A!(fi 2 ) 





M 



\ 








A 2 (0 2 ) 








r (-^) 1+7 1 






M(n 2 ) 


(-z) 1+27 








ml 


") 




ml 








'■■) 





The numbers Aj are the exact eigenvalues of M. Remarkably, they are determined by the 
part of dK/dx, which does not increase the degree in x, that is, by the operator on the 
left hand side of (|5TT|) : 



x n (n 2 ) 



7T 



-(O 2 + cos Tca n ). 



(5.3) 



sin na n 

where a n = 1 + rvy. By adequately adjusting O 2 , we can set any given A n to zero, which 
would correspond to the solution of the wave equation. Therefore, the spectrum is given 
by 

fi 2 = — COS7TCK n = COS7T717, (5-4) 
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so that 

P z ~ D — cos7rn7. 

If D < 1, then D = cos7T7, and P 2 ~ cos7T7 — cos7rn7. We see that, for 7 irrational, 
— cos7rn7 can get arbitrarily close to 1, thus again providing the evidence for tachyons in 
the spectrum. 

The only possible loophole in this argument could be the locality of our analysis. We 
did not find the global solution, but rather expanded it near the edge singularity. 

In fact, we can decide whether our solution corresponds to a minimum by checking 
that the eigenfrequencies of all normal modes are real and not imaginary for any value of 
O 2 . On the other hand, it is easy to see that the squares of eigenfrequencies of the operator 
—dK/dx have the same sign as A(0 2 ). If a n = 1 + 77-7, 

A n (Q 2 ) = (O 2 - cos7rn7). 

sin 71727 

Consequently, for small n these are always positive and do not cause any instability. How- 
ever, at large enough 717 > 1, when sin imy < the sign in front of P 2 changes, which 
creates ghosts. 

Although the higher n might indeed bring in ghosts and/or tachyons, one can imag- 
ine the situation where these modes are excluded by the boundary conditions. Strictly 
speaking, it is necessary to know the master field throughout the whole region of support 
to determine the spectrum fully. 



6. Conclusions. 

We have found that there is a special scaling domain where the shape of the master 
field in "induced QCD" is a universal function. We have demonstrated how to calculate this 
function in terms of a power series. Furthermore, we have constructed the D — 1-expansion 
to investigate the global structure of this power series. Finally, we have evaluated the meson 
spectrum, corresponding to this master field. 

Independently of the above, we have used the description of "induced QCD" in terms 
of the Hopf equation to derive a new, explicit, version of the wave equation for the meson 
spectrum of the theory. The equation we have found is valid for "induced QCD" with 
any potential and is not restricted to the particular application we have considered in this 
paper. 
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